Inclusive electron scattering off 4 He 



MKPH-T-07-02 



Sonia Bacca a *, Hartmuth Arenhovel fc , Mr Barnea c , Winfried Leidemann d and Giuseppina Orlandini d 
a Gesellschaft fur Schwerionenforschung, Planckstr. 1, D-64291 Darmstadt, Germany 
b Institut fur Kernphysik, Johannes Gutenberg- Universitat, J.J.-Becher-Weg 45, D-55099 Mainz, Germany 
c Racah Institute of Physics, Hebrew University, 91904, Jerusalem, Israel 

d Dipartimento di Fisica, Universita di Trento and INFN 
(Gruppo Collegato di Trento), via Sommarive 14, 1-38100 Trento, Italy 
(Dated: February 7, 2008) 

Inclusive electron scattering off 4 He is investigated for low and medium energy and momentum 
transfers. The final state interaction, given by the simple semirealistic Malfliet-Tjon potential, 
is treated rigorously applying the Lorentz Integral Transform (LIT) method. Besides the nonrel- 
ativistic one-body current a consistent meson exchange current is constructed and implemented. 
Results are presented for both longitudinal and transverse response functions at various momentum 
CD transfers. Good agreement with experimental data is found for the longitudinal response function, 

while some strength is missing in the transverse response function on the low-energy side of the 
quasi-elastic peak. 
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I. INTRODUCTION 



' In recent years considerable progress has been achieved in rigorous microscopic nuclear structure calculations of 
few-body nuclei with realistic or at least semirealistic interactions. Quite different approaches have been devised for 
O^l ■ bound state problems (see For break-up observables the powerful method of the "Lorentz Integral Transform" 
(LIT) has been developed in Ref . Q ■ The method is very convenient for calculating cross sections of perturbation 
induced reactions, where the inelastic response of the nuclear system is given by an absolute square of a transition 
matrix element summed over the whole spectrum of final states. Within this approach one does not evaluate the 
response function straightforwardly but calculates instead an integral transform of the response with a Lorentzian 
kernel. The response function itself is then obtained by inverting the transform. The specific advantage of this method 
lies in the fact that no explicit calculation of the complicated spectrum of final state wave functions is required. 

Over recent years, the LIT method has been exploited very successfully in the calculation of electromagnetic 
inclusive responses of various light nuclei (see, e.g. 0,11, H|), of neutrino inelastic reactions @,0|, and even of exclusive 
processes (see, e.g. With respect to inclusive electron scattering considered here, up to now only the longitudinal 

response function Rl, induced by the electro mag netic charge p(q), has been investigated in the LIT approach, namely 
for 3 H and 3 He with realistic nuclear forces [Tof, [ll| and for 4 He with the semirealistic TN potential [HI, [r|, while 
the transverse response Rt, induced by the current density operator j(q), has not yet been considered within this 
, method. 

In the pres ent work we evaluate both inclusive responses of the four-body system 4 Hc. With the exception of 
Rcfs. fl2l Il3ll . most previous calculations of Rl and Rt relied on the plane wave impulse approximation (PWIA) 
(see, e.g. |14|). where the interaction among the final reaction products is neglected (PW) and where only one-body 
electromagnetic operators arc considered (IA). In Ref. [HI, [l6[ the Laplace transforms of the two responses (the 
so called "Euclidean responses") were calculated using the Green Function Monte Carlo (GFMC) method with a 
realistic interaction and a consistent treatment of additional meson exchange currents (MEC) . This work showed that 
i) the final state interaction (FSI) is relevant in light nuclei, especially in the longitudinal response, and ii) two-body 
currents are important in the reaction mechanism, since they were found to have an enhancing effect of about 20% 
on the Euclidean transverse response. The strong FSI effect on Rt, of 4 He was confirmed in the above mentioned LIT 
calculation [l2l. [l3j with a semirealistic force, where reliable results were obtained not only for a transform of Rl, but 
also for Rt, itself. 

The main aim of the present work is the calculation of Rt for 4 He with a complete treatment of the final state 
interaction and using a current which is consistent with the potential. Since it is the first calculation using the 
LIT method we have chosen a semirealistic potential, i.e. the Malfliet-Tjon potential [17j . The longitudinal response 
function of 4 He will be reconsidered with this potential as well. In our calculation the final state interaction of the 
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four-body continuum is fully taken into account via the LIT method. Compared to the Laplace transform, the use of 
a Lorentz kernel has the advantage that a stable inversion of the integral transform is possible, thus allowing a direct 
comparison of the theoretical predictions with experimental data. 

In the following section we give a brief theoretical overview of the LIT method and remind the reader of the 
definition of the inclusive electron scattering response functions. Then in Sect. IIIII we discuss our results for the 
longitudinal and transverse response functions, and in the last section we draw some conclusions and give an outlook. 



II. THEORETICAL FRAMEWORK 



In this section we outline the main theoretical concepts for the present calculation of the electromagnetic response 
functions which govern the inclusive electron scattering off nuclei. First we give a short account of the Lorentz integral 
transform (LIT) which allows one to calculate the response functions with complete inclusion of the nuclear interaction 
in the final states. Then, in the second part, we briefly review the definition of the electromagnetic response functions 
and the model of the nuclear one- and two-body currents. 



A. The Lorentz integral transform 

As already mentioned, the basic observable is a response function which in the laboratory system has the general 
form 



V + — Ea-Lu) , 

} 2M ' ' 



where M is the target mass, I'J'o//) an d ^o// denote initial and final state wave functions and energies, respectively, 
while uj and q are the energy and momentum transfers. The operator O, describing the perturbative excitation 
mechanism, will be specified in the following subsection. The (^-function ensures energy conservation. The response 
function includes a sum over all possible final states, which are excited by the electromagnetic probe. This includes 
also the scattering states in the continuum. Therefore, in a straightforward evaluation one would need to calculate 
besides the initial ground state the final excited bound and continuum states of which the latter constitute the major 
obstacle for a many-body nucleus if one wants to include the nuclear interaction rigorously. 

In the LIT method [2| this difficulty is circumvented by considering instead of the response function Ro(uj,q) an 
integral transform Lo(<J, q) with a Lorcntzian kernel defined for a complex parameter a = an + iai by 

Lo(cr, q) = / dcu- . 2 

The imaginary part of a is kept at a constant finite value (07 7^ 0) determining the resolution of the integral transform. 
The basic idea of considering the integral transform is that it can be evaluated from the norm of a state \^„ q ), i-e. 

L (a,q) = {*° q \*° q ), (3) 
where \^„ q ) is the unique solution of the inhomogeneous "Schrodinger-like" equation 

(H-E -a)\*° q )=O(q)\y ). (4) 

Here H denotes the nuclear Hamiltonian. Because of the presence of a nonvanishing imaginary part 07 in Eq. ((D) and 
the fact that its right-hand side is localized, one has an asymptotic boundary condition like for a bound state. Thus, 
one can apply bound-state techniques for its solution. The response function itself, as a function of the ener gy t u for 
fixed q, is then obtained by inverting the integral transform for which various methods have been devised [l8l . . 



B. The response functions of electron scattering 



In the onc-photon-exchange approximation, the inclusive cross section for electron scattering off a nucleus is given 
in terms of two response functions, i.e. 



d 2 a 
dflduj 
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(5) 
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where <7m denotes the Mott cross section, Q 2 = —q^ = q 2 — lo 2 the squared four momentum transfer with u> and q 
as energy and three-momentum transfers, respectively, and 9 the scattering angle. The longitudinal and transverse 
response functions, Rl(w, q) and Rt{uj, q), are determined by the transition matrix elements of the Fourier transforms 
of the charge and the transverse current density operators, p{q) and j T (q), respectively, 




1. Multipole expansion 



It is useful to decompose the charge and current densities into Coulomb, longitudinal, and transverse electric and 
magnetic multipolcs |20j 



Jfj. v 



(8) 
(9) 



where {e M ; fi = 0, ±1} is a set of orthogonal spherical unit vectors with eo along q. Here the Coulomb and transverse 
multipole operators are defined as 



<tf(g) = dq'~p{q')Y^q'), 



I 

47T 

i 

47T 



dq'(q-3(q')K J (q') 



T raa 9 .J [q) = J d q< j {q >) . Y^(q>) 



d q '(q'xY^(q>)).j(q'), 



(10) 

(11) 
(12) 

(13) 



Coulomb and longitudinal multipolcs are related via current conservation, i.e. 



L J (q) = - H,C J Aq) 



(14) 



2. Siegert theorem and Siegert operators 

As is well known, the Siegert theorem states that in the low energy limit the transverse electric multipoles are related 
to the Coulomb multipoles. This is a very important and also very useful theorem, since it enables one to calculate, 
in the low energy regime, electric transition matrix elements from the charge density, without explicit knowledge of 
the current operator. 

Even beyond the low energy regime, one can always write each electric multipole as a sum of a Siegert operator 
plus a correction term. This is achieved by casting the transverse electric multipole operator in (Tj"2")) into the form 



2^(9) = ~J dq' 

= S^+K^q), (15) 
where the first term is related to the the longitudinal multipole operator (see Eq. (jll| ) 

S ^ q) = -^^T I dq'(q-j(q'W^q') = i]m^^( q ) (16) 



\f^-(q-3(q')K J (q') + 4jY J J+lfl (q') • j(q') 
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and 



ICfe) = --L-^= J dq'Y J J+1 ^q')-j(q'). 



(17) 



Using the relation (| 14[) from current conservation, the first term can be expressed by the Coulomb operator of Eq. IjlOl 
yielding 



J+i r 



H,C J Aq) 



(18) 



This form of S^(q) is called the li Siegert operator''. As is often said, the advantage of the Siegert operator is that, 
if evaluated with the one-body charge density, it includes already the dominant part of the MEC. The additional 
term K^(q) is a correction to the Siegert operator, and in the limit that the photon momentum goes to zero, the 
Siegert operator dominates the electric multipole, since the correction K'^(q) is two orders in q higher. Therefore, the 
approximation of the transverse electric multipolcs by the Siegert operators is quite reliable at low photon momentum, 
i.e. for qR <C 1, where R characterizes the spatial extension of the system. However, with increasing q, one has to 
calculate also the contribution of KAq), which requires the knowledge of the explicit form of the current operator. 



3. The current density operators 

In general, the electromagnetic current density of a nucleus can be decomposed in a superposition of one- and 
many-body operators. In a non-relativistic approach, as is the case here, the electromagnetic one-body four-current 
is given by the free nucleon current density operator, i.e. 

P{i){x)= e -Y,^ + rt)5(x-r k ), (19) 

k 

for the charge density, where is the third component of the isospin of the "fc-th" particle. The current density 
consists of a convection and a spin current 

= + (20) 

with 

f(i)(x) = ^^(l+T fc 3 ){p fe ,5(a;-r fc )}, 
k 

3h(x) = 5> s + r,V) <r k x \p k ,S(x - r k )} , (21) 

k 

where fi v , er k , Pki an d rn denote isoscalar and isovector nucleon magnetic moments, spin, momentum and mass of 
the "fc-th" particle, respectively. These expressions describe point particles. In order to take into account the internal 
nucleon structure, charge and magnetic moments have to be replaced by the corresponding form factors 

i(l + 4) - G|(Q 2 ) + 4GUQ 2 ) , M s + M V| - G S M (Q 2 ) + r fc 3 G^(Q 2 ) . (22) 

For on-shcll particles, these form factors depend on the squared four momentum transfer Q 2 alone. In principle, 
this is no longer true for the off-shell situation. However, in view of the fact that little is known about the off-shell 
continuation and furthermore for the moderate energy and momentum transfers considered here, the neglect of such 
effects is justified. Here we use for all form factors a common dependence on Q 2 of the usual dipole form [2lj | . 
In order to satisfy the continuity equation 

V-j i2) (x) = -i[V,p (1) (x)}. (23) 

a momentum and/or isospin dependent two-body interaction V requires a two-body current density operator, the 
interaction (or meson exchange) current. As is well known, relation (|23|) is not sufficient to determine the two-body 
current uniquely. Therefore one needs in principle a dynamic model for the nuclear potential which reveals the 
underlying interaction mechanism. Such a model is supplied by the meson exchange picture of the AW-interaction. 
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But even for a phcnomcnological potential often a consistent MEC can be constructed if the potential has the formal 
appearance of a meson exchange potential. This is achieved using the method developed independently in [22| and 
[23]. The method is based on the observation that if a potential is given as a sum of Yukawa- like terms, one can 
interpret the potential as produced by the exchange of fictitious mesons whose exchange currents are known. In the 
present case we use the semirealistic Malfliet-Tjon potential [l7j], which is isospin dependent and thus requires a MEC. 
Despite its phenomenological character, it fulfills the just mentioned requirement and thus it is possible to derive a 
consistent MEC which satisfies (J23J) with this potential. 
Explicitly, the MT-potcntial has the form 

V(r) = Via (r)Pi3 + Vai (r)P 3 i , (24) 

where Pi 3 and P31 denote the projectors on the spinsinglet-isotriplet and spintriplct-isosinglet two-nucleon states, 
respectively, 

P13 = -o-2)(3 + ti -Ta), (25) 

16 

P31 = 4( 3 + cr i • <T 2)(1 -Ti ■r 2 ). (26) 
16 

Furthermore, the radial functions arc given by 

V 13 {r) = i7r(AJ mi (r) - BJ m2 (r)), (27) 
Vai(r) = A7r(CJ mi (r)-DJ m2 (r)), (28) 

with J mi/2 (r) = e- mi / 2 747rr, ra a = 3.11 fin -1 , m 2 = 1.55 fin" 1 and A = C = 1458.047 MeV fm, B = 520.872 
MeV fm, D = 635.306 MeV fm. Exhibiting explicitly the spin-isospin dependence, one can rewrite the potential in 
the form 

V(r) = V (r) + V a {r) v 1 -tr i + (V T (r) + V aT (r) a x ■ <t 2 )t 1 ■ r 2 , (29) 

with 

V (r) = ^(^AJ mi (r)-BJ m2 (r)+CJ mi (r)-DJ m2 (r)y (30) 

V a (r) = ^(-3AJ mi (r)+WJ m2 (r)+CJ mi (r)-DJ rn2 (r)y (31) 

V T (r) - ^ (AJ mi (r) - BJ„ l2 (r) - 3CJ mi (r) + 3D J„ l2 (r)) , (32) 

Var{r) - ^(-AJ mi (r)+BJ m2 (r)-CJ mi (r) + DJ m2 (r)y (33) 

Because the isospin independent part commutes with the charge operator, only the isospin dependent part of (|29[) is 
relevant for the MEC. One can bring it into the form 

V T (r) + Vcrir) = 4ttti • r 2 (aJ mi (r) + (3J m2 {r) + (7J TOl (r) + 8J m2 (r))<Ti ■ cr 2 ) , (34) 

where the new constants are defined as 

a = i(A-3C), p = ~(B-3D), 
16 16 

1 = ~{A + C), 5=^{B + D). 

Now, since J m (r) represents the propagator of an exchanged scalar meson of mass m, the corresponding meson 
exchange current is easily constructed (for details see Ref. [2J]) and one finds 

J(2^' I ~ m ( a; ' r i) r 2) = 47r ( T i x x 2)3 aJ mi (\n - x\) V x J mi (|x -r 2 \) +0J m2 (\ri - x\) V x J ma (|x-r 2 |) 



-{lJ mi {\ri - x\) V x J mi (\x - 1-2 1) + SJ m2 (\n - x\) V x J m2 {\x - r 2 \)j &i ■ a 2 ■ (35) 
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It is a simple task to prove that this current fulfills the continuity equation (|23|) with the Malfliet-Tjon potential in 
For later purpose we list also the Fourier transform of the above MEC 



•MTI-III/ „ „ \ e 
3(2) (9^1> r 2) = 



iq-R 



-(ri x r 2 ) 3 aV r I mi (q,r) 



(3V r I m2 (q, r) + ax ■ ct 2 ( 7 V r / mi (<?, r) + 5V r I m2 (g, r) 



(36) 



where r and R denote relative and center-of-mass coordinates of the two-body sub-system, respectively, defined as 

r = n - r 2 and R=^(r 1 +r 2 ), 

and 



Jm 1/2 (9,r) = / d 3 p 



[(p + k) 2 + ™ 2 /2 H(p - h) 2 + TO i /2 ] 



(37) 



Here, we would like to recall that in contrast to a pseudoscalar meson exchange current the current in Eq. (|36|) does 
not include any contact term, but only meson in flight contributions, due to the fact that the MTI-III potential 
represents the exchange of scalar mesons. The functions Im 1/2 contain, in fact, the propagators of scalar mesons with 
effective masses mi and vn-i (for details see Ref. [H[). 



III. RESULTS AND DISCUSSION 



Now, we present results for the longitudinal and transverse electric and magnetic response functions for 4 He. We 
have solved the LIT-equation ([4]) by using an expansion of the internal wave functions in terms of hypcrspherical 
harmonics (HH). For the antisymmetrization we have used the powerful algorithm developed by Barnca and Novosel- 
sky [26[ . The HH-expansion has been truncated beyond a maximum value K max of the HH grand-angular momentum. 
Furthermore, the convergence of the HH expansion has been improved considerably by introducing the above men- 
tioned effective interaction with hypcrspherical harmonics (EIHH- approach) [27J. where the original bare potential has 
been replaced by an effective potential constructed via the Lee-Suzuki method [28|]. When convergence is reached, the 
results agree with the ones obtained with the bare potential [13 ■ As original potential we have used the Malfliet-Tjon 
MTI-III potential [13 • We found that for the ground state good convergence is reached already with K max = 10, 
yielding a binding energy of Eq = —30.57 MeV (the small difference with respect to the previously published value 
in [53 is due to the different K max used). For the state Vf^ °f Eq. (|3|) we have used K max = 11 or 12 depending on 
the odd or even parity of the excitation operator 0, respectively. The results we are presenting arc fully converged in 
K m ax- We begin the discussion with the longitudinal response function. 



A. The Longitudinal Response Function 



Since the longitudinal response function is an incoherent sum of the separate multipole contributions it is useful to 
expand the charge operator into Coulomb multipoles as defined in JHJ). For each multipole operator the LIT-equation 
with 0(q) = CjH(q) is solved. In view of the fact that isospin is a good quantum number and that the ground state 
of 4 He has total isospin zero, one can treat the isoscalar and isovector parts of the charge operator separately. The 
multipole expansion has been truncated at a maximum value J m ax determined by the requirement that convergence 
is achieved. It is known that with increasing momentum transfer J max has to be increased too in order to reach 
convergence. 

In Fig. [TJ we present the response functions of the isovector multipoles, obtained from the inversion of the LIT of 
each multipole, as a function of the intrinsic energy u>i n t = Ef — Eq. We show results for two momentum transfers 
q = 300 and 500 MeV/c. One readily notes that for the lower momentum transfer five multipoles are sufficient in 
order to reach convergence, while for the higher momentum transfer two additional multipoles need to be included. 
The dominant feature is the quasi-elastic peak located at uji n t — q 2 /2m, where m is the mass of the nucleon, with a 
shoulder on the low energy side. The latter decreases rapidly with increasing momentum transfer. 

The corresponding isoscalar response functions arc displayed in Fig. [5] The contributions of the isoscalar monopolc 
and of the sum of all multipoles up to J max are shown separately. An interesting feature of the monopole is the 
pronounced peak close to threshold which is even higher than the quasi-elastic peak for the lower momentum transfer, 
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FIG. 1: Response functions of the lowest isovector Coulomb multipoles, starting with the monopole and consecutively adding 
higher multipoles up to J ma x = 4 for q = 300 MeV/c (upper panel) and Jmax = 6 for q = 300 MeV/c (lower panel), as a 
function of the intrinsic energy. 

but suppressed for the higher momentum transfer, although still visible as a shoulder. In order to better determine the 
width of the monopole excitation of 4 He a more detailed study of low-g transitions with a realistic nuclear interaction 
is necessary. 

After the inclusion of the nucleon form factors and summing up the isoscalar and isovector contributions, we obtain 
the total longitudinal response function for varying energy transfers in the laboratory frame, where the recoil energy of 
the nucleus has been included. In Fig. [3j we compare our results with the available experimental data for momentum 
transfers q = 300 and 500 MeV/c. We have taken into account the Darwin- Foldy lowest order relativistic correction of 
the one-body charge density by a proper modification of the nucleon charge form factor [2^] . While the correction is 
negligible for q = 300 MeV/c, it leads to a damping of about 6% of the total longitudinal strength in the quasi-elastic 
peak for q = 500 MeV/c. The less important spin-orbit relativistic correction has been neglected (see e.g. fiojl)- 

As in the previous calculation of Rl using the Trento (TN) potential [lH, one notes that a semirealistic interaction, 
in this case the MTI-III, leads to quite a good overall description of the response in comparison to the experimental 
data from Bates [3(3 and Saclay [3l| . 

The only difference to the previous calculation with the TN potential is the pronounced peak close to threshold 
in case of q = 300 MeV/c which originates from the monopole excitation of 4 He. On the other hand, such a peak 
is not seen in the data. But it is not clear whether the experimental energy resolution was sufficient to resolve such 
a structure. It is worthwhile to mention that a + resonance at 20.10±0.05 MeV with a width of 270±50 keV was 
determined in an electron scattering experiment at momentum transfers q < 100 MeV/c (32j . Here we do not calculate 
these low-g kinematics, the resonance is very close to the "quasi-elastic peak" and quite small in size in comparison. 
A much more detailed study than the present calculation would be necessary to resolve such a rather complicated 
low-energy structure. 

B. The Transverse Response Function 

As done for the charge operator, we have expanded the transverse current operator into electric and magnetic 
multipoles according to Eqs. (fT2"| and (flU|) . separating them further into isoscalar and isovector parts, because the 
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FIG. 2: The response functions of the isoscalar Coulomb monopole and the summed multipoles up to J m ax = 4 for q — 
300 MeV/c (upper panel) and J max = 6 for q = 300 MeV/c (lower panel) as a function of the intrinsic energy. 

response function is an incoherent sum of these various multipole contributions. As discussed above, the transverse 
current includes one- and two-body operators. We first consider the one-body current alone, i.e. the spin and the 
convection current of Eq. (|2ip . Later we will add the consistent two-body current. 

1. One-body current 

It is known from standard PWIA calculations, that the spin current dominates the transverse response function 
at medium momentum transfers in the region of the quasi-elastic peak. Therefore, we start the discussion of the 
transverse response function of the spin current alone. In Fig. |4l we present the isoscalar and isovector response 
functions of the magnetic and electric multipoles up to J m ax = 4 for the spin current at a momentum transfer 
q = 300 MeV/c. One readily notes that the dominant isovector contribution exhibits a substructure in the quasi- 
clastic peak region, which originates mainly from the magnetic quadrupole and electric dipole contributions. A very 
pronounced peak is also found in the magnetic and electric isoscalar responses close to threshold, which is, however, 
suppressed in relation to the isovector part of the total response function after the inclusion of the nucleon form 
factors, because of the smallness of the isoscalar magnetic form factor relatively to the isovector one. 

The second contribution of the one-body current to the transverse response function is given by the convection 
current of Eq. (|2 1 1) . In this case one has a derivative term caused by the dependence on the nucleon momentum p k 
in the anticommutator, which acts on the wave function. In the EIHH approach this derivative can be calculated 
analytically, due to the fact that the radial wave function is given in terms of Laguerre and Jacobi polynomials with 
hyperradial and hyperspherical variables, respectively (see Ref. (27j). In order to prove the correct implementation 
of the derivative term, we make use of the Siegert operator as a check. Firstly, we write the convection current in 
momentum space as 



3(i) (?) = f(i)(l) +f(?)(Q), 



(38) 
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FIG. 3: (Color online) Longitudinal response function with the MTI-III potential as function of the ener gy f or momentum 
transfers q = 300 MeV/c (upper panel) and 500 MeV/c (lower panel). Experimental data from MIT-Bates [3[| (open circles) 
and Saclay [3lJ (filled circles). 



where 



4m 



qr k 



Pk 



(39) 
(40) 



In this way, the first part j c ,\\{q) does not contain any derivative term and is purely longitudinal, while the second part 



'(i) 

contains a longitudinal and a transverse, term depending on the direction of p k with respect to q. The check 
consists in the comparison of the Siegert operator of (| 18[) with the operator in (fT6)) using the convection current of (I38[) . 
For the case of a purely central potential, where no meson exchange currents are present, the two contributions have to 
coincide. As such an interaction we take for a test the model MTV, a modified version of the central Malfliet-Tjon-V, 
as defined in (24[. 

The results of this check for the dipole response are exhibited in Fig. [5] for different momentum transfers. First of 
all one can note that the dipole strength is differently distributed depending on the momentum transfer: for the two 
lower q- values the response function has a pronounced and quite small peak located a few MeV above threshold, while 
for the highest momentum transfer, q = 300 MeV/c, the peak becomes considerably broader, but smaller in height 
and shifted to higher energies. This different behavior has its origin in a different dependence of the two contributions, 
(?) an d 3°A\ (<?), to (fll)|) on the momentum transfer. The current (q) is dominant at q = 3 MeV/c, while for the 
two higher q- values also the second component, j'fft(q), becomes sizeable and interferes constructively with j c ^(q) at 
q = 100 MeV/c but destructively at q = 300 MeV/c. However, their total contribution agrees perfectly well with the 
Siegert dipole response. 

Again, the convection current is expanded into electric and magnetic multipoles, and then the response function 
is calculated for each multipole until convergence is reached. For example, in case of medium momentum transfers 
like q — 300 MeV/c we obtain a rather fast convergence pattern with J max = 4, since already only the dipole 
and quadrupole contributions cover about 80% of the total strength of the convection current. We also observe 
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FIG. 4: Transverse response functions at q = 300 MeV/c as functions of the intrinsic energy of various isovector (upper panels) 
and isoscalar (lower panels) magnetic (left) and electric (right) multipoles of the spin current, starting with the dipole and 
consecutively adding higher multipoles up to J ma x = 4. 



that only the electric multipoles are relevant, while the magnetic multipoles are negligible. For example, in case of 
q = 300 MeV/c, the magnetic dipole induced by the convection current is 10 times smaller than the electric dipole. 

Summing up the total contributions of spin and convection currents, we obtain the transverse response function in 
impulse approximation. In Fig. [6l we present the results for the MTI-III potential in comparison to the experimental 
data from Bates |3(| and Saclay [U] for momentum transfers q = 300 and 500 MeV/c. We also show the response 
given by the spin current only, which clearly dominates for both values of momentum transfer. For q = 300 MeV/c 
the addition of the convection current produces an enhancement of 3 — 4% near the peak energy uj = 60 McV. In 
case of q = 500 MeV/c, the convection current contribution is negligible. Qualitatively, the IA results reproduce the 
position and the width of the quasi-elastic peak, but show some deviations in detail. For example, for q — 300 MeV/c 
in the region of the maximum and above the theory is systematically below the data: in the maximum it is about 
10% lower compared to the Bates data and 20% with respect to the Saclay data. On the low energy side the theory 
exhibits a steeper rise than the data and in addition a secondary maximum which is not seen in the Saclay data, but 
might possibly be indicated by the Bates data. We would like to point out that evidence for an M2 resonance at 24 
MeV was found in the above mentioned (e, e')-experiment [32|]. On the other hand, for q = 500 MeV/c, the theoretical 
height of the quasi-elastic peak is only about 5% lower with respect to experiment, but the peak position is shifted 
by about 6 MeV towards higher energies. Furthermore, the theoretical shoulder near the threshold is barely seen in 
the data. These substructures at low energies arise from the electric dipole and magnetic quadrupole contributions 
of the spin current. 

It is conceivable that the overall missing strength could be provided by the up to now neglected two-body currents. 
Another explanation could be related to some inadequacies of the semirealistic description of the nuclear Hamiltonian 
or to rclativistic effects. In the following subsection we will investigate the first possibility. 
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FIG. 5: (Color online) Siegert part of the dipole response function for different values of the momentum transfer and evaluated 
with (i) the Siegert operator of Eq. (|18jl and (ii) the operator of Eq. (|16p with the convection currents j'^(q) and j%.(q), and 
their sum (see text) for the MTV potential. 



2. Two-body current 

The consistent two-body current for the MTI-III potential, j^J l ~ lU (q, r i, r%) is given in (PH)) . Its general multipole 
decomposition is listed in Eq. (|A10[) of the appendix. In the present MEC calculation we neglect the center of mass 
motion of the two-body subsystem, i.e. we set e lq R = 1. We would like to point out that in this approximation all 
magnetic and even electric multipoles vanish identically (see appendix). One expects this approximation to be quite 
good for small momentum transfers, i.e. qR ^ 1. In Fig. we compare, for the most important dipole contribution, 
results with explicit MEC and convection current on the one hand and with the Siegert operator on the other hand. 
In contrast to the MTV case depicted in Fig. \5\ with the MTI-III potential, the Siegert dipole does not agree with 
the response of the Siegert part of the convection current alone, due to the presence of the MEC. At low momentum 
transfer the convection current leads to a response, whose peak height is only about half of the Siegert peak, while 
at higher momentum transfer the effect of MEC becomes even stronger. The dipole response function of the MEC 
is also shown in Fig. [JJ It is interesting to note that the strength of the convection current is located near the 
quasi-elastic region whereas the MEC has a much broader distribution. Taking as excitation operator the sum of 
convection current and MEC, we get a perfect agreement with the Siegert dipole for q = 3 MeV/c, which clearly 
proves that the approximation introduced in the two-body current is reliable at low momentum transfer, as expected. 
For q = 100 MeV/c our approximation can still be considered very good. However at the higher values of momentum 
transfer, one notes an increasing difference between the two evaluations which has to be assigned to the neglected 
cm. motion of the two-body subsystems in the calculation of the MEC. One further observes that the effect of the 
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FIG. 6: (Color online) Transverse response function in impulse approximation (IA) and response of the one-body spin current 
alone as a function of the energy for momentum transfer q — 300 (upper panel) and q — 500 MeV/c (lower panel) in comparison 
with data from Bates [33] and Saclay [3ll |. 
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FIG. 7: (Color online) Dipole response function with the MTI-III potential, as a function of the energy and for different values 
of the momentum transfer. The response is evaluated for (i) the Siegert operator of (|18[) (solid curves) and for the operator of 
(|16|) with (ii) convection current (dotted), (iii) MEC (dash-dotted) and (iii) their sum (dashed). 



approximation consists mainly in an enhancement the peak height of the response, while position and form of the 
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peak are not affected. The peak reduction amounts to about 10% and 20% for q = 200 and 300 MeV/c, respectiveiy. 
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FIG. 8: Dipole response function with the MTI-III potential as a function of the energy, for different values of the momentum 
transfer. The response is evaluated for (i) the Siegert operator of (|18l) (solid) and (ii) for total dipole operator according to 
(|15|) with the sum of convection current with MEC (dashed). 



Furthermore, we can check to what extent the Siegert operator alone is a good estimate of the total dipole operator, 
i.e. whether the correction term of (fT7|) can safely be neglected. To this purpose we consider the total multipole of 
Eq. (fj"5j) . where the correction term (fT7|) has been included in addition to the Siegert part (fT6|) . and we compare it to 
the Siegert operator itself alone. In Fig. [51 we show the result of this check again for the isovector dipole response as 
a function of the energy for different momentum transfer values. The total dipole response function induced by the 
convection current and by the MEC with mutual interference leads to the same result as the Siegert dipole response 
for the very low momentum transfer q = 3 MeV/c, and even up to a momentum transfer of q = 100 MeV, the Siegert 
operator is still a reliable approximation. This is not anymore the case at higher momentum transfers where one 
readily notes an increasing difference between the total dipole response and the Siegert one alone. For momentum 
transfers q = 200 and 300 MeV/c, the Siegert operator clearly underestimates the total dipole response at energies 
below the peak and overestimates it in the tail. In fact, the strength of the Siegert response is moved as a whole 
towards higher energies with a shift of the peak position with respect to the total dipole response by about 17 MeV 
and 33 MeV for q = 200 and 300 MeV/c, respectively. On the other hand, with respect to the total dipole response 
the peak height of the Siegert response is reduced by about 10% at q = 200 MeV/c and enhanced by about 20% for 
q = 300 MeV/c. Obviously, the Siegert operator is not a good approximation for the total dipole response for q > 100 
MeV/c. At ui = 100 MeV, for example, compared to the result of the total dipole response, the Siegert response leads 
to an overestimation of 35% and 70% for q = 200 and 300 MeV/c, respectively. 

Though the situation may be quite different for other electric multipoles, like quadrupole and octupole, and even 
more for the isoscalar multipoles, we think that the results of this check on the dominant isovector dipole clearly 
suggests that one needs to consider the total electric multipoles when investigating the transverse response function 
at medium to high momentum transfers, since in this kinematic region the Siegert approximation is not reliable 
anymore. Thus, an explicit knowledge of the two-body current operator is necessary. 

In Fig. El wc finally present the transverse response function for different momentum transfers in comparison to the 
available experimental data. We compare the IA with the full calculation in which we add the two-body currents. For 
q = 100 and 200 MeV/c we also show the response function induced by the spin current only. In Table U we recall 
explicitly the multipoles considered in our final calculations of Rt{lo, q) for the different transverse current parts and 
for each momentum transfer value. As previously shown in Fig. [HI the spin current strongly dominates the response 
function at a momentum transfer q = 300 MeV/c. A rather similar situation is found for q = 200 MeV, as depicted in 
Fig. [9j On the other hand, one can see that for the lower momentum transfer q = 100 MeV/c the convection current 



14 




FIG. 9: (Color online) Transverse response function for different momentum transfer values: spin current only, IA and full 
calculation with the inclusion of a consistent two-body current, in comparison with the available experimental data from 
Bates [3(| and Saclay [3l|. (For q = 300 MeV/c the result for the spin current only is already shown in Fig. [BJ. 



TABLE I: Maximal multipolarity J m ax considered in the calculation shown in Fig. [9] for the one- and two-body currents and 
for the different momentum transfer values. 
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presents a stronger effect: for example, at ui = 40 MeV it is even half of the total IA strength. By comparing the IA 
with the full calculation, we can conclude that quite a strong effect of the MEC is found at low momentum transfers 
like for q = 100 MeV, where for example at energies of 100 MeV they lead to an enhancement of the strength by a 
factor of four with respect to the IA. At this momentum transfer, the approximation we had introduced for the MEC 
is safely reliable. As the momentum transfer increases, for example at q — 200 MeV/c, the enhancing effect of the 
MEC with respect to the IA result on Rt is reduced dramatically and it amounts to only 4% in the peak region and 
to 40% for u> = 100 MeV. Unfortunately, to the best of our knowledge, no experimental data have been measured in 
the quasi elastic region at these low momentum transfers q = 100 and 200 MeV/c. For the case of q = 300 MeV/c 
where experimental data from Bates [3(| and Saclay [3l[ exist, we observe a tiny effect of the two-body current with 
respect to IA result. The addition of the MEC leads to an increase of 4% and almost 10% at w = 100 and 150 MeV, 
respectively. Our full calculation agrees rather well with the data of Bates in the energy range 55 < ui < 115 MeV 
but is lower than the data of Saclay in the quasi-elastic peak region and above. For example, at ui = 60 MeV our full 
result for q = 300 MeV/c is still lower than the measurement of Saclay. 

Here we would like to mention that the approximation introduced in the MEC by neglecting the center of mass 
dependence, which is estimated to be of about 10% to 20% in the peak and negligible in the tail for q = 200 to 
300 MeV/c, does not affect our conclusion, due to the fact that the overall MEC effect is small. 
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IV. CONCLUSIONS AND OUTLOOK 



We have presented the first calculation of the inclusive longitudinal and transverse response functions of 4 He within 
the LIT and EIHH methods with complete inclusion of the interaction in the final states as given by the semirealistic 
potential MTI-III potential. As in previous LIT calculations with the semirealistic TN potential [12| a good overall 
agreement with available experimental data is found for the longitudinal response function, though the MTI-III 
potential leads in addition to the quasi-elastic peak to another peak close to threshold, which is not clearly seen in 
the data of Bates [13] and Saclay [HI, while, as mentioned, at q < 100 MeV/c a resonance close to threshold has been 
observed in experiment [32| . 
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FIG. 10: (Color online) Euclidean transverse response for momentum transfer q — 300 MeV/c. Comparison of the GFMC 
calculation [l(| (band between thin lines) and the result of this work (thick line) with the Laplace transform of experimental 
data from Bates [3(| and Saclay [3l| . 



For the transverse response function, we have performed a calculation in IA alone and then adding in a second 
step a two-body current consistent with the MTI-III potential. Strong MEC effects are found at low momentum 
transfer, q — 100 MeV/c, where unfortunately no experimental data are available. For the case of q — 300 MeV/c, for 
which data from Bates [3(| and Saclay [3l| exist, we have shown that the IA result still misses some strength in the 
quasi-elastic peak region. For the present MEC, consistent to the semirealistic potential model, we have not found a 
strong two-body current effect at q = 300 MeV/c as obtained in the calculation of Carlson and Schiavilla within the 
Laplace transform approach [l6[ including it- and p-MEC. A direct comparison of the two calculations for R,T{uj,q) 
is not possible because the inversion of the Laplace transform suffers from large instabilities. In Fig. 1101 we thus 
compare the Laplace transform of our full transverse response function for q = 300 MeV/c with the GFMC Euclidean 
resp onse, using as interaction the Argonnc V18 potential [H| with the addition of the Urbana IX three-body force 
[34| . We adopt the definition of the Euclidean response as presented in [ll| [l6[ and we integrate our theoretical 



curve up to the maximal value of energy transfer (u) max = 180 MeV) measured in Saclay [3l| for the considered 
momentum transfer case. Consistently with what is obtained in a direct comparison, our calculation for the MTI-III 
potential with the present MEC model leads to a lower transverse Euclidean response with respect to the Laplace 
transform of experimental data. The Euclidean response obtained with the GFMC, where contributions of 7r- and 
p-exchange currents were considered, agrees better with data. We have to conclude that the additional contribution 
of the present MEC does not explain the missing strength of our model. The difference to the GFMC calculation 
could possibly be due to the missing explicit pionic degrees of freedom in our potential model and therefore also in 
the corresponding MEC. In particular, our MEC model does not include any contact terms because the fictitious 
meson exchanges comprise solely scalar mesons. It is well known, that such contact terms as associated with pion 
exchange currents are quite important. Moreover, the masses of the exchanged mesons of the present potential model 
are considerably higher than the pion mass making the exchange current shorter ranged and thus less important in 
the range of momentum transfers considered in this work. Thus it might not be surprising that we find a reduced 
effect of MEC. On the other hand it is interesting to notice that the IA result for MTI-III is closer to the data than 
the corresponding one with AV18+UIX. This lets open the possibility that for MTI-III the relative contribution of 
two-body currents is indeed smaller than for a realistic interaction. Future studies should be extended to realistic 
potential models including pion degrees of freedom. 
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APPENDIX: MULTIPOLE EXPANSION OF THE MEC 



Here, we will give a brief derivation of the multipole expansion of a two-body current of the form 

3(2){q,r,R) = e^ R j (2) (q,r) 



(Al) 



consisting of an intrinsic part j ( 2 ), which depends only on the relative two-body coordinate r, and a two-body cm. 
part, depending on the two-body cm. coordinate R. The consistent meson exchange current of ([36|) has exactly this 
form. 

The basic quantities to evaluate are 



J J L,mM = ^ / dq'e^ R Y J LM (q') ■ Vj (2) (q\r) , 



(A2) 



where Y J L M denotes a vector spherical harmonics. The quantities in (|A2I) determine the longitudinal and transverse 
electric and magnetic multipolcs according to 
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el. J 



L 

L=J±1 
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L 1 J 
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J L.O i 



LI J \ t j 
1-1 Lt " 



J 1 J 
1 -1 



Jj _ jj 



(A3) 
(A4) 
(A5) 



for = 1. Here a spherical coordinate system with fi £ {0,±1} has been chosen where eo is along q. For more 
details the reader should consult Ref. [24| . 

In order to evaluate (|A2[) one hrst expands e lq R and the intrinsic current J(2)(q, v) into spherical harmonics. The 
former one is given by the well known expression 



e iqR = 4ir^2(-) jR J R \c Jn (q, R) x Y Ja (q) 
Jr 

where we have introduced the cm. Coulomb multipolcs 

C J *{q,R)=i J *jj R {qR)Y jR {R), 
with spherical Bcssel functions ji(x). The corresponding expansion for the intrinsic current is given by 



Jj.inj.Lj. 



where 



J 



'Jr 

Lj.,m 



in 



dq'Yt, mr (q')-V'j {2) (q',r) 



(A6) 



(A7) 



(A8) 



(A9) 



Inserting these expansions into (|A2[> . then, after recoupling, the angular integration can be done analytically, and one 
finds finally 



JiJq,r,R) = (-) J+1 V^JL ^ J R L r 

J r J r L r 



Jr L r L \ J Jji J r J 







1 L L r 



C J «(q,R)xjl : (q,r) 



(A10) 
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With the help of angular momentum algebra, the longitudinal and transverse multipoles are then obtained from (|A3 
through (|X5]) . i.e. 



Li( g ,r,R) = (-y^Lj2^ 



J 



JuJr 



Jr J r J 




C J *(q,R)xL J *(q,r) 



(All) 



•J r„ j \ / 



J R Jr 

C jR (q,R) X f el / ma9 ' J r(q 7 r] 



(A12) 



Thus a total longitudinal or transverse multipole is given by a sum of all possible couplings of corresponding intrinsic 
and cm. charge multipoles as allowed by parity and angular momentum coupling rules to form a rank- J tensor. 

With respect to the multipole expansion of the consistent MEC of ([55)) we start from the expansion V/ m , which is 
given in [2 51 ] 



with 



VI m (q, r) = 4tt J2 J £m(9> r ' m ) Y tM 

J r fiL r —even 



1 L J 



If^q, r, m) = 4 7 r(i) J - 1 L f J J q ) Y ^ r ' m) ' 

where the functions $^ ; (g,r, m) arc defined by 
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1 JO 



with z — (p 2 + \q 2 + m 2 )/pq and 



Qt{z) = - 



-jcr(pr)Qi(z), 



dx- 

1 z — X 



(A13) 



(A14) 



(A15) 



(A16) 



Here, Pi{z) denotes a Legendre polynomial. Then one obtains for the intrinsic part of the MEC in 

^J r ,L r (<7: r : m i) + P < S'j r ,L r ( ( l> r > rn 2) +<T i ' ^(t^^l,.^,^"*!) +S<S> Jr Lr (q,r,m 2 



, (A17) 

which, inserted into (|A11[) and (|A12|) finally yields the MEC multipoles. From the condition L r = even follows in 
conjunction with the 3j-symbol in (|A14[) . that j£ r (q,r) vanishes for L r = J r and consequently all intrinsic magnetic 
and even electric multipoles vanish according to (|A4[) and (|A5|) . 

If one neglects the cm. contributions, i.e. setting R = 0, then one has 



JL(Q,r,0) = JUq,r). 



(A18) 



This approximation is used in the present work. 
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